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We investigate the threshold of gravitational collapse with angular momentum, under the assump- 
tion that the critical solution is spherical and self-similar and has two growing modes, namely one 
spherical mode and one axial dipole mode (threefold degenerate). This assumption holds for perfect 
fluid matter with the equation of state p — up if the constant k is in the range < k < 1/9. There is 
a region in the space of initial data where the mass and angular momentum of the black hole created 
in the collapse are given in terms of the initial data by two universal critical exponents and two 
universal functions of one argument. These expressions are similar to those for the correlation length 
and the magnetization in a ferromagnet near its critical point, as a function of the temperature and 
the external magnetic field. We discuss qualitative features of the scaling functions, and hence of 
critical collapse with high angular momentum. 
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I. INTRODUCTION 



An isolated system in general relativity ends up in one of three final states. It either collapses to a black hole, 
forms a star, or disperses completely. The phase space of isolated gravitating systems is therefore divided into basins 
of attraction. One cannot usually tell into which basin of attraction a given data set belongs by any other method 
than evolving it in time to see what its final state is. The study of these boundaries in phase space, in particular of 
the boundary between black hole formation and dispersion, is the subject of the new field of critical collapse ||l|,^. 

The pioneering work of Choptuik |^ has shown that the black hole threshold is both richer in structure and simpler 
than naively expected. Choptuik carried out systematic high precision numerical collapse simulations in the toy model 
of a spherically symmetric massless scalar field coupled to general relativity. He explored the black hole threshold 
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by means of smooth one-parameter families of initial data. A generic family contains both collapsing and dispersing 
data sets. Choptuik found that locally each family crosses the black hole threshold only once. This suggests that the 
collapse threshold is a smooth hypersurface in phase space. 

The black hole mass shows a universal power-law scaling as a function of distance from the black hole threshold, 
as measured by the parameter p of the one-parameter family of initial data. Let a black hole be formed for p > p^, 
and let dispersion occur for p < p^. The critical value p^ of the parameter p depends on the family, and can only be 
determined through a bisection search. Choptuik found that the black hole mass is approximately 

M{p) ^C{p-p,y' forp>p„ (1) 

where the critical exponent 7 is independent of the family (universal), with a numerical value of 7 ~ 0.374 for the 
scalar field. 

Furthermore, all evolutions with p ~ p^, on either side of p*, and for any family, pass through an intermediate 
attractor, before finally dispersing or forming a black hole. This "critical solution" is self-similar. In the case of the 
scalar field, the self-similarity is a discrete symmetry. In the case of the perfect fluid, it is continuous. By virtue of 
being an intermediate attractor, the critical solution has precisely one growing mode. The critical exponent for the 
black hole mass can be calculated from this growing mode. 

Subsequent work has found similar critical phenomena in many other matter models (see for a review). Unfor- 
tunately, most work so far has been limited to a spherically symmetric, uncharged situation in which a Schwarzschild 
black hole is formed. Generic black holes, however, have angular momentum L and electric charge Q as well as mass 
M, and we do not yet know what happens if one fine-tunes to the black hole threshold along a one-parameter family 
of data with significant angular momentum and/or charge. 

Today, the role of angular momentum and charge in critical collapse is understood only in the limit where they are so 
small that they can be treated as linear perturbations of a spherically symmetric and uncharged scenario throughout 
the collapse, from the initial data to the final black hole. One can then make a connection from a small perturbation 
in the initial data representing charge or angular momentum to a linear perturbation of the final Schwarzschild black 
hole that takes it into a Kerr-Newman black hole. 

Specifically, critical collapse with a small amount of electric charge has been investigated in the model of a spherically 
symmetric complex scalar field coupled to a Maxwell field. A critical exponent for the black hole charge Q has been 
calculated by keeping track of the least slowly decaying charged perturbation of the critical solution Q . This prediction 
has subsequently been verified in collapse simulations Q . 

Similarly, critical collapse with a small amount of angular momentum has been investigated in the model of a 
perfect fluid with equation of state p = Kp, where p is the pressure, p the total energy density, and k is a constant 
with < K < 1. In exact spherical symmetry, the critical exponent 7 for the black hole mass is again independent 
of the initial data but depends on k [^-^ . Perturbing around spherical symmetry, a critical exponent for black hole 
angular momentum L was calculated from the least slowly decaying axial dipole perturbation in [ic[| . The numerical 
value of the critical exponent has been corrected in . For the scalar field, a critical exponent for L was calculated 
in second order perturbation theory around spherical symmetry ]T^ . Both these predictions still await testing, as 
high-precision simulations of rotating collapse to a black hole are not yet available. 

In this paper we go beyond the assumption that angular momentum is a small perturbation throughout. This is 
necessary because the corrected results for the axial dipole perturbations of the critical fiuid in show that all such 
perturbations decay only for the equations of state k > 1/9. For k < 1/9, there is actually a single growing mode. 
This means that in the limit in which initial data with a small amount of angular momentum are fine-tuned to the 
black hole threshold, there are two competing growing modes of the critical solution, one spherical and one to do 
with differential rotation. The latter is threefold degenerate because rotation is an axial vector. Which of the two 
modes first reaches nonlinearity depends on their initial amplitude and growth rate. In particular, if the initial data 
are sufficiently close to the collapse threshold, rotation may dominate the late stages of the evolution even if it was 
only a small perturbation in the initial data, and the final black hole could be very small but rapidly rotating. 

We begin our presentation in section || by re view ing the qualitative picture of critical collapse in terms of a dynamical 



system, and establishing notation. In section III we identify an intermediate attractor through which solutions near 



the black hole threshold are funnelled. In section tV we discuss the initial data that reach this attractor, and in 
we follow them from the intermediate attractor to the black hole end state. We find that the black hole mass 
and angular momentum depend on the initial data only through universal functions of one combination of the initial 
data parameters. In this paper, we do not calculate these "universal scaling functions", but their mere existence 
substantially constrains the phenomenolo gy of critical collapse resulting in small but rapidly spinning black holes. 



We discuss this phenomenology in section VII. The present work has been motivated by the close analogy between 



critical collapse and critical phase transitions in sta tistical mechanics, and re pre sents an extension of that analogy. 



This analogy is discussed separately in section VIII. We summarize in section IX. 
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In talking about black holes, it is customary to refer not to its angular momentum L but to its specific angular 
momentum a = L/M, and we adopt this convention here. 



II. THE DYNAMICAL SYSTEMS PICTURE 



The time evolution of general relativity can be considered as an (infinite-dimensional) dynamical system once the 
coordinate freedom of general relativity has been fixed by imposing suitable coordinate conditions on parts of the 
metric. Here we review this phase space picture in preparation for the calculation of critical exponents and scaling 
functions in the next section. Black holes on the one side, and flat spacetime on the other, are attracting fixed points. 
Their basins of attraction are separated by the black hole threshold, or critical surface. It is obviously a hypersurface 
of codimension one. The numerical evidence is consistent with the assumption that it is smooth, and not, for example, 
fractal. By definition, the critical surface is a closed dynamical system in its own right. Within the full phase space, it 
is a repeller. Its attracting fixed points or limit cycles are then attractors of codimension one in the full phase space, 
with exactly one growing perturbation mode. They are called critical fixed points in dynamical systems language, or 
critical solutions in spacetime language. 

Any trajectory beginning near the critical surface, but not necessarily near the critical point, moves almost parallel 
to the critical surface towards the critical point. As the critical point is approached, the parallel movement slows 
down, and the phase point spends some time near the critical point. Then the phase space point moves away from the 
critical point in the direction of the growing mode, and finally ends up on a stable fixed point. This is the origin of 
universality: any initial data set that is close to the black hole threshold (on either side) evolves to a spacetime that 
approximates the critical spacetime for some time. When it finally approaches either empty space or a black hole it 
does so on a trajectory that appears to be coming from the critical point itself. 

Not surprisingly, attractors of this dynamical system have higher symmetry, when considered as Cauchy data for a 
spacetime, than generic points in the phase space. The critical solutions found at the black hole threshold, at least 
those known to date, are either static, periodic, continuously self-similar, or discretely self-similar. Here we consider 
only the continuously self-similar critical solutions, which give rise to scaling laws of the form (|l|). The model we focus 
on here, the perfect fluid with equation of state p = np, has this type of critical solution. A spacetime is continuously 
self-similar, or homothetic, if there exists a vector field x such that the Lie derivative of the spacetime metric along 
X obeys 

i^xQab = -Sfifah- (2) 

X is called a homothetic vector field. It is a special case of a conformal Killing vector field. By the Einstein equations, 
continuous self-similarity implies that the matter stress-energy obeys C^Tab — 0. 

In general relativity, each data set locally determines a unique spacetime. But a given spacetime can be broken 
up into a sequence of data sets, that is, foliated by spacelike hypersurfaces, in an infinite number of ways, even if 
the initial hypersurface is kept fixed. This means that the same physical spacetime can be described by an infinite 
number of curves through the phase space of general relativity. Turning general relativity into a dynamical system 
therefore requires some prescription that returns a lapse and shift for each point in phase space. In the following, let 
r be the time coordinate of the dynamical system, and let x be the three spatial coordinates. 

Calculations in critical phenomena require a prescription in which the self-similar spacetimes are fixed points, and 
discretely self-similar spacetimes are limit cycles. This means that for those sets of Cauchy data that will evolve into 
a self-similar spacetime, the prescription returns a lapse a and shift vector /3" such that the time vector d/dr is the 
homothetic vector x- 

|:) ^an^ + (3^^x'', (3) 
where n° is the unit normal vector on the constant r slices. Certain possible prescriptions, such as "i^-freezing lapse" 



with "minimal distortion shift" have been identified |13|, and we now assume that one such choice has been made. 
Note that while d/dr becomes spacelike at large distances from the center of collapse, this does not preclude surfaces 
of constant r from being everywhere spacelike. The spacetime metric in these coordinates is 

9tiu{T,x) Z^e"^^5^j,(x), (4) 

where g^j/ does not depend on r, and where I is an arbitrary fixed length scale. Foliating by surfaces of constant r 
we obtain an induced metric and intrinsic curvature of the form 
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gij{T,x)=fe '^''gijix), Kij{T,x) = le ''Kij{x). (5) 



We now turn this around to define the rescaled dynamical variables 

gij{T,x) = r'^e^'' gij{T,x,), Kij{T,x) = r'^e^'' Kij{t,x) (6) 

for the metric, and similar rescaled variables for the matter fields. For a perfect fluid with 4- velocity and comoving 
energy density p, for example, we define 

u' = le-''u\ p = fe-^''p. (7) 

In the variables Z = {gij, Kij, m*, p, . . .}, and with a suitable choice of lapse and shift, a CSS solution will take the 
form Z(x, t) — Z^(x). Full physical initial data consist of Z{x, tq) and the value of tq itself, which provides an overall 
scale and allows one to reconstruct the physical variables from the barred variables. 

It is helpful to think of the coordinates x and t and the barred variables as dimensionless, of Ze^^ as having a 
dimension of length (or time, or mass, in units c = G = 1), and of the unbarred quantities as having their natural 
dimensions, with g^^ having a dimension of because the coordinates are dimensionless. The Einstein and matter 
equations in the barred variables are much like the usual ones, but are dimensionless. e^'^ only appears in the 
combination le~'^ . In vacuum gravity, with a minimally coupled massless scalar field, or with a perfect fluid with 
equation of state p = up fluid, no dimensionful constants appear in the field equations. Therefore there are no 
constants to cancel Z, and so no powers of e^"^ can appear explicitly in the dimensionless equations in the barred 
variables. This means that solutions are possible which do not depend on r, and so are self-similar. If dimensionful 
constants of dimension Z"", with n > 0, appear in the field equations (for example, a mass term in the scalar wave 
equation), they are multiplied by Z"e^"'^ in the scale-invariant equations, and so become dynamically irrelevant as 
T —f oo, which corresponds to very small scales. This allows for solutions that become asymptotically self-similar on 
small scales (large t). 



III. THE INTERMEDIATE LINEAR REGIME 

For what follows it is helpful to reformulate the mass-scaling law in spherical symmetry, Eq. (|l|), by absorbing the 
family-dependent constants C and p* into a new parameter p that is a linear function of the parameter p, so that the 
law is now 

M for p > (8) 

for every 1-parameter family of initial data. It is an experimental fact that this is possible, with each data set assigned 
only one value of p, independently of the 1-parameter family of which it is considered a part. This means that p is a 
scalar function on the infinite-dimensional phase space, independently of any 1-parameter families. M is also a scalar, 
but it is not regular at the black hole threshold, while p is regular with non-vanishing gradient, p is therefore a good 
coordinate on phase space in a neighborhood of the black hole threshold. This clarifies the invariant meaning of the 
critical exponent 7: it is the unique power that relates the observed scalar M on phase space to another scalar p that 
is a good coordinate. 

Based on the phase space picture, the critical exponent 7 for the black hole mass can be calculated in a manner 
suggested by Evans and Coleman and spelled out by Koike, Hara and Adachi [Q and Maison Here we repeat 
this calculation, but now including the effects of angular momentum in the initial data. We do this in the abstract 
notation Z which applies not only to the perfect fiuid with k < 1/9, but equally to any other model that has a 
spherically symmetric, continuously self-similar solution with precisely two growing perturbation modes. One of these 
must be spherically symmetric, and the other one must have an axial dipole (1 = 1) angular dependence. All other 
spherical and nonspherical perturbation modes (there are infinitely many) must decay. The generalization from a 
continuously self-similar to a discretely self-similar critical solution is trivial. 

As the background solution we consider is continuously self-similar, with Z{x,t) = Z^,{x), the perturbations of 
this background depend on r only exponentially. We can write a generic perturbation on the background Z^{x) as 
a sum over perturbation modes of the form e^''^Zi{x, il), where the index i labels the perturbation spectrum. The 
two growing modes in particular are designated Zo{x) (for the spherically symmetric one) and Zi{x,il,) (for the axial 
dipole mode). Here the arrow indicates that Zi is 3-fold degenerate, for m = — 1,0, 1. fl is shorthand for the angular 
dependence on the coordinates {0,ip). Perturbing around spherical symmetry, to linear order rotation is associated 
with the axial dipole mode, and in particular, the angular momentum of the spacetime is proportional to the 
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fallofF of the dipole mode at spacelike infinity. This is true for regular perfect fluid spacetimes, but also for black 
holes: outside the horizon, a Kerr black hole with a <^ M can be written as an axial dipole linear perturbation of 
a Schwarzschild black hole. The significance of the axial dipole perturbation is discussed in more detail in Refs. [T^ ] 
and | [l5| ]. As a matter of notation, we go from the basis m = —1, 0, 1 to the basis xyz of / = 1 harmonics, so that Zi 
is an axial vector in 3-space that points in the direction of angular velocity. 

In this paper, we consider solutions that go through an intermediate time regime in which they are well approx- 
imated by a self-similar, spherically symmetric solution plus small perturbations, the "intermediate linear regime" . 
Furthermore, in this regime, we neglect the decaying linear perturbation modes, so that we have the approximation 

Z{x,T,Vl) ~ Z^{x) + Ae^"^ Zfi{x) + B ■ e^^^ Zi{x,9) + Aec&yi^g modes. (9) 
In the following we show that all solutions in a ball straddling the hole threshold are in fact funneled through the 



intermediate linear regime. Working backwards, in section IV we construct families of non-spherically symmetric 
data that go through the intermediate linear regime. Working forwards, in section we extract Cauchy data in 
the intermediate linear regime, and predict the mass and specific angular momentum of the black hole from these 
intermediate data. 



IV. INITIAL DATA WITH ROTATION 

We now construct generic families of initial data that reach the intermediate linear regime (|). To do this, we shall 
first introduce a parameter controlling angular momentum, and then add a second, generic, parameter which allows 
us to tune to the black hole threshold. 

Consider first infinitesimally small deviations from a spherically symmetric data set that does form a black hole. 
(For the moment we do not assume that we are at the black hole threshold.) Perturbations with different multipoles, 
and axial and polar perturbations, then remain decoupled throughout the evolution, and so the angular momentum 
of the final black hole is linearly related to the axial dipole part of the initial data. However, while the specific angular 
momentum of the black hole is a vector a (three numbers) , the axial dipole part of the initial data for a perfect fluid 
consists of a vector of functions /3(r), which characterize differential rotation. We can explore this 3oo-dimensional 
space by means of many different 3-parameter families of initial data, with the vector-valued parameter q chosen so 
that q — Q corresponds to spherical symmetry, and so that q —q corresponds to /3(r) —P{r) for all r. Clearly, 
this is a sufficient (but by no means a necessary) criterium for the relation a,{—q) = —a{q) in the final black hole. It 
is also clear that we have M{—q) = M{q) for the black hole mass. (In fact, from this it follows that to linear order in 
perturbation theory, M does not depend on q. 

But we need not stop at perturbations of spherical symmetry. In order to make Schwarzschild black holes, one does 
not need spherically symmetric initial data. Note that a Schwarzschild black hole has three commuting reflection 
symmetries (also called octant symmetry), while a Kerr black hole has only one reflection symmetry (through the 
equatorial plane). Any reflection symmetry in the initial data is maintained during evolution. Octant symmetry 
in the initial data is therefore a sufficient (but not a necessary) criterium for producing a non-rotating black hole. 
This leads us to consider 3-parameter families of initial data that are not a perturbation of spherical symmetry: we 
construct them so that the data with q = Q are octant-symmetric, and that data with q^O contain a point such that 
q —q corresponds to a reflection through that point. This symmetry between q and —qis conserved by the time 
evolution, and so again we have the relations 

a{-q) = -a{q), M{~q) = M{q). (10) 

While octant symmetry of the initial data is a sufficient criterium for the absence of angular momentum, it seems 
unlikely that one can give a necessary criterium for the final black hole to have zero angular momentum. But we can 
turn the argument around and formally define the parameter q to have the properties (p^. Such a definition may 
seem circular, but it is in fact the analogue of the definition of the scalar function p, in that from the family-dependent 
parameter q we can define a parameter q that is a scalar on the space of initial data and which is related in a simple 
universal way to the black hole specific angular momentum a. 

We have defined the vector parameter q without assuming that we were close to the black hole threshold. However, 
if the black hole threshold is a surface of codimension one in the full space of non-spherical initial data, then fine- 
tuning any other parameter p in the initial data should take us to the black hole threshold, for any fixed value of q. 
With this motivation, consider now 4-parameter families of initial data with the following properties: 

1. The initial data depend analytically on p and g in a neighborhood of the black hole threshold. 

2. If for any p and q a black hole is formed, its specific angular momentum and mass obey ( |lO| ) for any p. 
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3. For (f= 0, there is a p* such that a black hole forms if and only if p > p*. 

We should note that is the value of p at the threshold only for (f = 0. The black hole threshold is locally described 
by P = Pciit(g) with Pcrit(O) = p*. 

Do data sets from such a family go through the intermediate linear regime (|9|)? 

a) li q — corresponds to spherical symmetry, and if l^l is sufficiently small, the answer is yes. In this case, the 
deviations from spherical symmetry are described by perturbation theory throughout. We can then appeal to the 
experimental fact that in spherical symmetry, the basin of attraction of the spherically symmetric critical solution is 
the entire black hole threshold, including data that are far from the critical solution. 

b) If q does not correspond to spherical symmetry, the answer is unknown. For example, there is no guarantee that 
the spherically symmetric critical solution is reached from near-critical data that are octant-symmetric (and therefore 
non-rotating) but highly oblate or prolate. But it is suggested by the fact that a spherical critical solution exists, and 
that critical solutions apparently have as much symmetry as possible. (In critical collapse in pure gravity, the critical 
solution cannot be spherically symmetric, but it is axisymmetric) . 



V. UNIVERSAL SCALING FUNCTIONS 

We are now ready to calculate a(p, g) and M{p,q). The amplitudes of Zq and Zi during the intermediate linear 
regime depend on p and q in some complicated, nonlinear, and family-dependent way. Assuming that the dependence 
is analytic, we can determine the behavior to leading order from symmetry considerations. From the definition 
and because Zi carries the angular momentum during the intermediate linear regime, B is an odd function of q, and 
A is an even function. Also, both A and B must vanish for p = and g = because by definition that initial data 
set is on the black hole threshold. To leading order in p and q, we find that 

Z{x,t) ~ + Co(p-p,)e^"^Zo(:E) + {Ciq) ■ e^'^ Zi{x,n). (11) 

Here Co and Ci are two unknown constants that depend on the two-parameter family we consider. Ci is a 3 x 3 
matrix. In axisymmetry it reduces to a number, but in general it is a non-trivial map between the vector q and the 
vector Zi, and therefore, in the end, between q and a. 
We now define r» such that 

Co|p-p,|e^«^* =e, (12) 
where e is a fixed small constant. We then have 

Z{x,T^) ~ Z^{x) ± eZo(x) +5- Zi{x,n), (13) 

where 



Here the sign in front of e is that of p — p^ . It appears because of the absolute value taken in the definition (|l^) . 
Recall that in the following e is a fixed positive constant, while 5 depends on the initial data. 

In order to simplify our notation, it is convenient to introduce "reduced" values p and q of p and q as 

p = Coe"^(p-p,), q = Ciq, (15) 

Note that by definition the direction of q is now the direction of a. In the following, this is implied, and we suppress 
the vector indices. In this simplified notation, 

S=\p\-^q, (16) 

The parameters p and q can be thought of as "reduced" values of p and q, which hide the family-dependent part of 
the dependence of M and a on p and q. They are similar to the "reduced temperature" and "reduced external field" 
in statistical mechanics. We can also think of them as coordinates on the phase space, in the way already discussed 
for p. 

We now consider Z(a:,T») as Cauchy data for a time evolution in r that leads away from the perturbative regime 
and eventually to black hole formation or dispersion. Z(a;,r,) contains the complete Cauchy data, up to an overall 
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scale, which must be provided separately. That scale is e""^*, which from our definitions is proportional to 

In particular, the mass of the black hole that is formed must be proportional to this scale, with a constant of 

proportionality that can depend only on S and the sign of p. We find 

M(p,,) = H*{g<*j; f/^l (17) 

The functions F^f{S) are universal, and so are the two exponents Ao and Ai. Because of ([l0|), Fm is an even function. 

Consider the special case q = 0, and hence 6 = 0. Then only the mode Zq is present in the intermediate linear 
regime, and, up to scale, we only need to consider the two data sets — eZq and + eZo- It is known that the first 
data set disperses, while the second one forms a black hole. If we use M = to denote the absence of a black hole, we 
have F^j{0) = and F^j-{0) > 0. We use the freedom to normalize Zo{x) in order to set ^"j^(0) = 1 as a convention. 
We then find 

M{p,q)^{p^^ P>0 asg-->0. (18) 

We have recovered the power law in its form (||), with the critical exponent 7 = I/Aq for the black hole mass. 

We now consider the specific angular momentum a of the black hole. The dimensionless quantity a/M can depend 



on the initial data only through the dimensionless quantity S, and on the sign of p. M itself is given by (17), and we 
only have to put these two results together to obtain a. It is convenient to absorb the scaling function Fm{S) into 
the new scaling function for a, and we can write 

.(P.,)-iPi-{?f<l!; (w) 

Here F^{S) are two new universal scaling functions. Clearly, -F'^(O) = 0. Recall that we have adapted a simplified 
notation, in which Fa and S are really vectors, but in which Fa takes its direction trivially from S, so that Fa is really 
just a scalar function of a scalar variable. Because of (p^), it is an odd function. If cosmic censorship holds, then 
a < M must hold in black holes formed in gravitational collapse, and therefore \Fa\ < Fm for all 5. 

To leading order Fa is proportional to 5. By adjusting the normalization of Zi, we can set the factor of propor- 
tionality to one. With this convention, and from the definition (p^, we obtain a power law for a in the limit of small 
T- 

( _ 

a{p,q)^\(iP ^ P>^ asg^O. (20) 
lO, p<0 

This had been obtained previously under the assumption that Ai < [|lO[. 



VI. PROPOSED NUMERICAL TESTS 



So far our results have been formal: we know the two exponents Ao and Ai, but not the universal scaling functions. 
However, we predict that such functions exist and are universal. This can be tested by comparing two or more two- 
parameter families of initial data. The two family-dependent parameters C'o /e and Ci (the latter is a matrix) must 
be determined for each family, by fitting observations in the limit g — > to the formulas (|l^) and (|2^). The scaling 
functions can then be read off from the first family, and tested against the second and any further family. 

Alternatively, the scaling functions can be calculated directly by evolving the data (^3|) for fixed e and all values of 
S. In the notation of [ll]] , the area radius r in spherical symmetry is related to the dimensionless radial coordinate x 
and the scale coordinate and coordinate r by r = sxe^"^ . (s is a known constant whose significance does not matter 
here.) We set r = 0, and consider the initial data 

^initiai(r,r!) = Z4r/s)±AZo{r/s)+BZi{r/s,n). (21) 

Here Z^{x) is already known, and Zo{x) and Zi(x) are known up to normalization ||ll|. From the three-fold degenerate 
Zi, we arbitrarily choose to = 0. The actual dimensionless variables that Z stands for in perfect fluid collapse are 
defined in pl] |. The physical, dimensionful fluid variables of are obtained from the rescaled variables defined in 
by setting r = everywhere. 
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We now evolve these data nonlinearly in axisymmetry. In order to explore all values of 5 in the universal scaling 
functions, we could keep A fixed at a small value ±e, and vary B from to oo. In practice, however, we want to vary 
the ratio A/B from —oo to oo, but choose an overall factor in both A and B small enough so that the initial data 
are in the intermediate linear regime but large enough so that the evolution leaves this regime soon. This minimizes 
the range of scales that the nonlinear evolution code has to cover, and may make the calculation possible without the 
need for adaptive mesh refinement. 

Clearly, p is proportional to A and q is proportional to B. The constants of proportionality, however, are unknown 
because we do not have the correct normalizations of Zq and Zi. We just fix an arbitrary normalization, and put in 
two adjustable constants. Then the mass and specific angular momentum of the final black hole that is created from 



the initial data (21) are related to the scaling functions by 

M{±aA,f3B) = A'^F^j (^BA^^^ 



(22) 



a{±aA,l3B)^A^F^ {BA ^) . (23) 



The constants a and /3 correspond to the unknown correct normalizations of Za and Zi, and must be adjusted to 
obtain the two conventions 

^^+(0) = 1, F+'(0) = 1. (24) 



VII. LARGE ANGULAR MOMENTUM 



Current axisymmetric rotating fluid evolution codes are not ready at the time of writing to calculate the universal 
scaling functions in the manner suggested in the previous section. It is therefore tempting to speculate about the 
form of the universal scaling functions at large 6. 

Consider the limit where p —>■ while g 7^ so that 5 ^ 00. In physical terms this corresponds to first fine-tuning a 
one-parameter family of initial data without angular momentum to the black hole threshold and then adding a small 
but finite amount of angular momentum to the data. Do the data obtained in this way form a black hole? If our 
model is correct, the answer is universally yes or no, independently of the family of zero angular momentum data we 
have fine-tuned and independently of how we have added that bit of angular momentum. The reason is that such 
data are funneled through the data set Z^ + eZ\ (for some small fixed e) in the intermediate linear regime. We now 
discuss the two possible final states for these data in turn. 

Possibility 1: The data -I- eZ\ form a black hole. As its mass is finite (these data possess an intrinsic scale). 



the explicit power of p in formula (17) must be canceled by a power-law behavior of Fj^i{S) as i5 — > 00. One easily 
sees that this power law must be 

F±(<5)~Cm|'5F, \S\^^ (25) 
for the cancellation to occur. But then one immediately obtains that 

M-CmI^F, \5\^(x. (26) 



By the same argument we have 



and therefore 



Ft (6) ^Ca signi6)\d\-, |<5|^oo, (27) 



a^CaSign{q)\q\-, \S\ ^ 00. (28) 

The constants Ca and Cm are again universal. We see that a/M goes to the constant Ca/CM as we increase the 
angular momentum in the initial data. We do not know what the value of this constant is, except that cosmic 
censorship requires Cq/Cm < 1- 

Possibility 2: The data -I- eZi do not form a black hole. On physical grounds this alternative appears more 
likely, as one would expect centrifugal forces to disrupt data that already hover between collapse and dispersion. 
This means that F^^(±oo) = 0. If we consider obtaining this limit by adding angular momentum to supercritical 
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data, that is by increasing q at fixed p > 0, it appears likely that above a threshold amount of angular momentum 
no black hole is formed at all, rather than ever smaller ones. In this case it also seems plausible that zero angular 
momentum data below the black hole threshold will not form a black hole if one adds angular momentum. If these 
two assumptions are true, then instead of (|25| ) and we have 

^M,a(^) = for \5\ > ,5,nax, Fj^j JS) = for all S. (29) 
This means that black holes are formed if and only if 

_ _ _— 

p>0 and |g| < (5inax ■ (30) 

This defines a convex region in the p, q plane. We do not know if black holes at the threshold with non-zero angular 
momentum are formed with zero or finite mass. This gives rise to two sub-cases: 

Possibility 2a: In the presence of angular momentum there is a mass gap at the black hole threshold. This means 
that F^j is discontinuous at ^max, of the form 

for some universal constant Km and the size of the mass gap is 

/^M~KmP'^> (32) 
Similar behavior would hold for , for some universal constant Ka so that 

a^a-Kap'^ (33) 
Therefore, the spin/mass ratio of black holes formed at the threshold is universal, with 

= (34) 
AM Km 

From cosmic censorship we have Ka/KM < 1. Fig. |^ gives a qualitative impression of M{p, q) and a{p, q) in case 2a. 

Possibility 2b: There is no mass gap even in the presence of angular momentum. This means that F^{S) vanishes 
at (5max- If this happens, for example, as a power n of distance from the threshold, 

then the black hole mass would scale as that same power n of any regular coordinate scalar on the phase space that 
vanishes at the threshold (except at q = 0). Fig. gives a qualitative impression of M{p, q) and a(p, q) in case 2b. 



VIII. THE ANALOGY WITH STATISTICAL MECHANICS 



This work was motivated by the attempt to exploit the critical phase transition/critical collapse analogy to learn 
something new about critical collapse. The preceding material has been presented in a self-contained manner without 
any explicit reference to statistical physics concepts, but it may be interesting to point out the exact parallels now. 

The calculation of critical exponents in critical collapse is mathematically identical to the calculation of critical 
exponents in statistical mechanics, even though the underlying physical phenomena are totally different. The mathe- 
matical equivalent of the renormalization group acting on micro-states in statistical mechanics is the time evolution, 
in certain preferred coordinates, acting on initial data in general relativity. Both can be considered as dynamical 
systems. Critical exponents are calculated by linearizing around their critical fixed points |p|,p^. 

In statistical mechanics, the critical fixed point typically has two growing modes. Physically these are linked to 
the temperature T and to a generalized external force. It is helpful to consider two contrasting examples. For a 
fluid confined in a vessel so that a vapor phase is in equilibrium with a liquid phase, the generalized force is the 
pressure P, and the order parameter is the difference puquid — Pgas between the densities of the two phases. For a 
ferromagnet, the generalized force is the external magnetic field H , and the order parameter is the magnetization 
rh. In both cases the temperature-force plane contains a line of first-order phase transitions ending in a second-order 
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phase transitions. (In the standard terminology, a first-order transition is one where the order parameter jumps from 
zero to a non-zero value. A second-order, or critical, phase transition is one where the order parameter is continuous, 
rising as a power-law.) For the fluid, this is the liquid-gas phase transition, ending at the critical point where these 
two phases have the same density. For the ferromagnet, the phase transition is between the possible directions of 
the magnetization m at T < Tc- This first-order transition ends at the Curie temperature Tc- The spontaneous 
magnetization at zero external field behaves as a power of Tc — T for T ^ Tc- 

If one restricts consideration to a vanishing external magnetic field, one appears to have an analog of the black-hole 
threshold in spherical symmetry, where the Tc — T corresponds to p — and the absolute value of the magnetization 
|m| corresponds to the black hole mass M. The formula corresponds to the dependence of the spontaneous 
magnetization on the temperature, in the absence of an external field. But typical critical phase transitions in 
statistical mechanics have two independent critical exponents, corresponding to two growing modes of the fixed point. 
Nigel Goldenfeld has asked if a second growing mode could not be found in the critical collapse problem so that the 
universal critical behavior would show not one critical exponent but two and, more interestingly, universal scaling 
functions 0- 

Here we suggest that the analog of the external magnetic field H in the collapse problem is the angular momentum 
parameter q. Angular momentum provides a second growing mode of the critical point. The analog of the order 
parameter to, a vector, is not the black hole mass but its (specific) angular momentum vector a. Conversely, the 
analog of the black hole mass M is a scalar that scales like a length, for example the correlation length ^. (The fact 
that ^ diverges at the critical point, while M goes to zero because that is what length scales do in the two kinds 
of critical point.) We have shown that a second critical exponent and universal scaling functions do indeed arise, in 
exact mathematical analogy with the ferromagnet model [ p^ . 

In both the ferromagnet and collapse models the two growing modes have different symmetries. While + eZg 
forms a black hole, Z^—eZg disperses. The two signs are qualitatively different. But Z^ + eZq + S- Zi will either collapse 
or disperse depending only on the absolute value of 6, independently of its direction. Similarly in the ferromagnet, 
the state at T = Tc + e differs qualitatively from that aX T = Tc — e, but only the absolute value of H matters. In 
each case, the vector-valued parameter {q or H) is associated with a symmetry breaking. Because q ~ or H = 
correspond to an unbroken symmetry, their critical value is obviously zero, whereas the critical values Tc and are 
nontrivial. {Tc depends on the material, and on the family.) If one imposes the symmetry on the problem (no 
external magnetic field, or octant symmetry in collapse), the critical point has effectively only one unstable mode. 

There is one major difference between the ferromagnet and critical collapse: In the ferromagnet, an infinitesimal 
external field at a temperature slightly above the Curie temperature creates a finite magnetization. The ferromagnetic 
region is therefore concave in the T-H plane. Even with zero external field, the symmetry is spontaneously broken, 
and a finite net magnetization in a random direction results. Fig. ^ gives a qualitative impression of m{Tc — T, H) 
for a ferromagnet in axisymmetry, and should be contrasted with Figs. |l| and A partial analog of this in critical 
collapse would be possibility 1 in the previous section. However, possibility 2 appears more likely on the physical 
grounds that centrifugal forces should oppose collapse. The black hole region is then convex in the pq plane. 

The analogy between critical collapse and the critical point of a fluid is less close than between critical collapse 
and the ferromagnet, as neither temperature nor pressure are associated with a symmetry breaking: both have to be 
fine-tuned at once to nontrivial critical values that depend on the material. As a curiosity, however, we note that 
an analog of the liquid-gas transition, where g is a scalar parameter, and not associated with a symmetry breaking, 
could arise in critical collapse in the context of semiclassical gravity. Self-similarity in critical collapse in semiclassical 
gravity is broken at small scales (of the order of the Planck scale) by the quantum stress-energy. Brady and Ottewill 
||18| suggest that the effect of the quantum stress-energy can be described approximately in terms of a single unstable 
mode (in addition to the known classical one) with a growth exponent of Ai = 2. If this description is correct, the 
parameter q would describe some measure of the excitation state of the quantum fields. 

IX. CONCLUSIONS 

The black hole threshold in gravitational collapse has been investigated until now only for initial data without 
angular momentum. It is largely unknown what happens when the threshold is crossed along a family of initial data 
with nonzero angular momentum. If the angular momentum is very small, it can be treated as a perturbation of 
spherical collapse throughout, and for this case a critical exponent for the black hole angular momentum has been 
and is now waiting to be tested in numerical collapse simulations. 
Here we have noted that the perturbative calculations giving rise to this prediction can be extended in the case where 
the angular momentum perturbation constitutes a second growing mode of the self-similar, spherically symmetric 
critical solution, in addition to the spherical growing mode known already. We can now extend our predictions to 
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a regime where the angular momentum of the initial data is small, but large enough so that the final black hole is 
rapidly rotating. 

In Eqs. ( p^ and (|l^) we have made predictions for the black hole mass and specific angular momentum in this 
regime, expressed in terms of two known critical exponents and two yet unknown functions of one variable S. The 
prediction that universal scaling functions exist at all can be tested by comparing two or more two-parameter families 
of initial data. The scaling functions could also be calculated more simply and precisely through the nonlinear 
evolution of only two one-parameter families of axisymmetric initial. 

The necessary initial data are already known to high precision, and the necessary technology (an axisymmetric fiuid 
code that can form rotating black holes, but without mesh refinement) should be available soon. In the meantime we 
have discussed the qualitative features of these functions, distinguishing two possibilities, and suggesting one of them 
as the most plausible on physical grounds. 

The prediction of universal scaling functions in critical collapse with angular momentum both extends and clarifies 
the analogy between critical phenomena in gravitational collapse and critical phase transitions in statistical mechanics, 
in particular the ferromagnetic phase transition. 
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FIG. 1. Schematic plot of black hole mass M and specific angular momentum a in axisymmetry as functions of p and q, 
assuming case 2a. For illustration we have assumed Ao = 2, Ai = 1, F+(5) = e{l - \5\) and F+(5) = 0(1 - \5\) sin(7r<5/2). Note 
that in reality Ai/Ao <IC 1, and that we have assumed that \a\/M^ — » 1 at the black hole threshold. 
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FIG. 2. Schematic plot of black hole mass M and specific angular momentum a in axisymmetry as functions of p and q, as- 
suming case 2b. For illustration we have assumed Ao = 2, Ai = 1, F+{S) = 61(1 - cos(7r5/2) a.nd F+{5) = 6»(1- sin(7r5)/2. 
Note that in reality Ai/Ao <C 1, and that we have assumed that \a\/M^ — > 1 at the black hole threshold. 
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1 ^ 

FIG. 3. Schematic plot of the magnetization m of a ferromagnet as a function of Tc — T and H. For illustration we have 
assumed Ao = 2, Ai = 1, F' {5) = sign(5)|<5|i/^i and F+{5) = F' {5) + e"^' [sign(5) + 8\. 
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